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1.  Introduction 

The  boundary  layer  of  a  mature  hurricane  has  been 
long  recognized  as  an  important  feature  of  the  storm.  In 
particular,  it  controls  the  radial  distribution  of  moisture, 
vertical  motion  and  absolute  angular  momentum  that 
ascends  into  the  eyewall  clouds.  Early  attempts  to  isolate 
the  dynamics  of  the  boundary  layer  of  a  mature  hurricane 
were  made  by  Rosenthal  (1962),  Smith  (1968),  Carrier 
(1971)  and  Eliassen  (1971).  These  studies  highlighted 
a  feature  already  well  known  from  other  areas  of  fluid 
dynamics  that  the  boundary  layer  induces  a  secondary 
(overturning)  circulation  in  the  vortex  above  and  that  this 
circulation  is  associated  with  the  imbalance  of  forces 
in  the  layer  brought  about  by  surface  friction  (e.g. 
Greenspan,  1968).  A  scale  analysis  of  the  boundary-layer 
equations  shows  that  the  vertical  gradient  of  perturbation 
pressure  can  be  neglected  to  a  first  approximation  so 
that  the  radial  pressure  gradient  in  the  boundary  layer 
is  equal  to  that  in  the  flow  above  the  layer  (Smith,  1968). 
However,  surface  friction  reduces  the  tangential  wind 
speed  and  hence  the  centrifugal  and  Coriolis  forces  in  the 
layer,  leaving  for  a  cyclonic  vortex  a  net  inward  residual 
force.  It  is  this  net  force  that  drives  inflow  in  the  layer, 
thereby  inducing  a  secondary  circulation  in  the  vortex 
itself.  Consistent  with  mass  continuity,  there  is  induced 
subsidence  at  outer  radii  and  induced  ascent  at  inner  radii. 
In  a  hurricane,  the  inflowing  air  acquires  moisture  through 
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evaporation  at  the  sea  surface  and  the  moist  air  ascends 
to  feed  the  eyewall  clouds. 

Despite  that  fact  that  the  induced  flow  in  the  bound¬ 
ary  layer  is  associated  with  gradient  wind  imbalance  in 
the  layer,  many  representations  of  the  hurricane  bound¬ 
ary  layer  have  taken  the  tangential  flow  there  to  be  in 
strict  gradient  wind  balance,  but  have  included  a  sink  of 
absolute  angular  momentum  at  the  surface  (Ogura,  1964; 
Ooyama,  1969;  Schubert  and  Hack,  1983;  Emanuel,  1986, 
1989,  1995,  1997,  2004;  Frisius,  2005,  2006;  Wirth  and 
Dunkerton,  2006).  We  are  unaware  of  any  reasoned  jus¬ 
tification  for  such  a  formulation,  but  one  could  imagine 
that  the  dynamics  of  the  balanced  boundary  layer  is  essen¬ 
tially  different  from  that  in  which  friction  plays  an  active 
role,  since  any  inflow  in  these  balance  models  cannot  be 
‘driven’  directly  by  frictional  imbalance.  At  best  they  may 
be  expected  to  be  a  valid  first  approximation  only  where 
the  radial  flow  is  sufficiently  slow.  A  natural  questions 
then  arises,  how  severe  are  their  limitations?  We  have 
shown  recently  (Smith  et  al. ,  2008)  that  the  assumption 
of  balance  is  poor  in  Emanuel’s  steady-state  hurricane 
model  (Emanuel,  1986)  and  by  implication  in  his  the¬ 
ory  for  potential  intensity  (Emanuel,  1995;  Bister  and 
Emanuel,  1998). 

Two  of  the  early  pioneering  models  of  hurricanes, 
those  of  Ooyama  (1969)  and  Sundqvist  (1970),  were 
balanced  in  the  sense  that  the  tangential,  or  primary, 
circulation  was  taken  to  be  in  axisymmetric  gradient 
wind  balance,  and  both  models  were  considered  to 
produce  reasonably  realistic  simulations  of  hurricanes. 
Indeed,  the  axisymmetric  gradient  balance  assumption 
is  thought  to  be  a  relatively  accurate  one  over  most  of 
the  free  troposphere,  except  in  the  upper-level  outflow 
layer,  and  the  assumption  is  supported  by  an  elementary 
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scale  analysis  for  a  rapidly  rotating  vortex  in  which 
the  radial  component  of  flow  is  much  less  than  the 
tangential  component  (Willoughby,  1979).  It  is  supported 
also  by  aircraft  reconnaisance  measurements  in  the  lower 
troposphere  (Willoughby,  1977,  1990).  However,  these 
early  models  assumed  that  the  boundary  layer  is  in 
gradient  wind  balance  also.  Ooyama  was  aware  of  the 
limitations  of  the  latter  assumption  and  wrote  in  an 
unpublished  manuscript  in  1968  ‘...  it  appears  that  the 
weakest  hypothesis  in  the  original  model  is  the  use  of 
the  balance  approximation  in  the  boundary  layer’.  In 
the  manuscript,  Ooyama  went  on  to  show  that  solutions 
with  a  more  complete  boundary-layer  formulation  were 
more  realistic  than  those  with  a  balanced  boundary-layer 
formulation.  As  far  as  we  are  aware,  most  authors  have 
not  seriously  questioned  the  accuracy  of  the  balanced 
assumption  in  the  hurricane  boundary  layer,  even  though 
some  prior  analyses  of  the  layer  (Smith,  1968;  Carrier, 
1971)  showed  that  it  is  unlikely  to  be  an  accurate 
approximation  in  the  inner  core,  a  finding  supported 
also  by  more  recent  numerical  simulations  of  hurricanes 
(Persing  and  Montgomery,  2003,  Appendix;  Kepert  and 
Wang,  2001). 

In  this  paper  we  revisit  this  important  problem.  We 
present  first  a  novel  scale  analysis  that  helps  identify  a 
baseline  of  rational  approximations  to  the  steady  depth- 
averaged  swirling  bondary  layer.  Then  we  compare  the 
predictions  of  various  approximate  formulations  of  the 
boundary  layer  in  a  steady  slab  model,  including  those 
that  assume  gradient  wind  balance,  with  that  of  the 
unapproximated  formulation.  Using  the  unapproximated 
solution  as  a  benchmark,  the  accuracy  of  the  various 
approximate  formulations  is  assessed. 

We  begin  by  reviewing  the  features  of  balanced  models 
in  general,  highlighting  the  role  of  friction  in  the  sec¬ 
ondary  circulation  as  described  by  the  Sawyer-Eliassen 
equation. 


2.  Balanced  hurricane  models 
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where  r,  k,  z  are  the  radial,  azimuthal  and  vertical  coordi¬ 
nates,  respectively,  (u,v,w)  is  the  velocity  vector  in  this 
coordinate  system,  C  =  v2/r  +  fv  is  the  sum  of  the  cen¬ 
trifugal  and  Coriolis  terms,  /  is  the  Coriolis  parameter, 
P  =  p/p  is  the  pressure  p  divided  by  the  mean  density  p 
at  height  z,  b  is  the  buoyancy  force  per  unit  mass,  defined 
as  —  g{p  —  p(z)}/p*9  p  is  the  density,  p*  is  the  average 
density  over  the  whole  domain,  B  is  the  diabatic  source 
of  buoyancy,  and  F,  and  F-A  are  the  radial  and  tangential 
components  of  frictional  stress,  respectively. 

With  the  additional  assumption  of  gradient  wind  bal¬ 
ance,  Equation  (1)  reduces  to 


(6) 


If  P  is  eliminated  by  cross-differentation  with  the  hydro¬ 
static  Equation  (3),  we  obtain  the  thermal  wind  equation 
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where  £  =  2v/r  +  /  is  twice  the  absolute  angular  veloc¬ 
ity  at  radius  r.  The  SE  equation  is  obtained  by  differ¬ 
entiating  (7)  with  respect  to  time,  eliminating  the  time 
derivatives  of  v  and  b  using  (2)  and  (4)  and  introducing 
a  streamfunction  xfr  for  the  secondary  circulation  such 
that  the  continuity  equation  (5)  is  satisfied,  i.e.  we  write 
u  =  —  (l/r)(dx/f  /dz)  and  w  =  (l/r)(dx/f  /dr).  Then,  with 
a  little  algebra  we  obtain 
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The  cornerstone  of  all  balance  theories  for  vortex 
evolution  is  the  Sawyer-Eliassen  (SE)  balance  model 
which  describes  the  slow  evolution  of  an  axisymmet- 
ric  vortex  forced  by  heat  and  (azimuthal)  momentum 
sources.  Here  the  flow  is  assumed  to  be  axisymmetric 
and  in  strict  gradient  wind  and  hydrostatic  balance.  We 
summarize  the  SE  model  here  for  the  simplest  configura¬ 
tion,  namely  the  axisymmetric  flow  of  an  incompressible 
Boussinesq  fluid  with  constant  ambient  Brunt-Vaisala 
frequency,  N.  The  hydrostatic  primitive  equations  of 
motion  may  be  expressed  in  cylindrical  polar  coordinates 
(r,  k,  z)  as 


where  S  =  dv/dz  is  the  vertical  shear  of  the  tangential 
wind  component  and  £a  =  (l/r)(drv/dr)  +  /  is  the 
absolute  vertical  vorticity.  More  general  derivations  of 
this  equation  are  found,  for  example,  in  Willoughby 
(1979),  Shapiro  and  Willoughby  (1982)  and  Smith  et 
al.  (2005).  The  SE  equation  is  elliptic  if  the  vortex 
is  symmetrically  stable  (i.e.  if  the  inertial  stability  on 
isentropic  surfaces  is  greater  than  zero).  It  is  readily 
shown  that  symmetric  stability  is  assured  when 
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(Shapiro  and  Montgomery,  1993).  Given  suitable  bound- 
^  '  ary  conditions,  Equation  (8)  may  be  solved  for  the  stream- 
function  of  the  secondary  circulation  xj/  at  a  given  time. 
(2)  Being  a  balance  model,  only  one  prognostic  equation  is 


Copyright  ©  2008  Royal  Meteorological  Society 


Q.  J.  R.  Meteorol.  Soc.  134:  1385-1395  (2008) 
DOI:  10.1002/qj 


BOUNDARY  LAYERS  IN  HURRICANE  MODELS 


1387 


used  to  advance  the  system  forward  in  time.  The  set  of 
equations  (2)  (or  (4)),  (7)  and  (8)  thus  provide  a  sys¬ 
tem  that  can  be  solved  for  the  balanced  evolution  of  the 
vortex.  Equation  (2)  (or  (4))  along  with  the  diagnostic 
equation  (7)  is  used  to  predict  the  future  state  of  the  pri¬ 
mary  circulation  with  values  of  u  and  w  at  a  given  time 
being  computed  from  the  streamfunction  \j/  obtained  by 
solving  (8).  The  secondary  circulation  given  by  (8)  is  just 
that  required  to  keep  the  primary  circulation  in  hydrostatic 
and  gradient  wind  balance  in  the  presence  of  the  pro¬ 
cesses  trying  to  drive  it  out  of  balance.  These  processes 
are  represented  by  the  radial  gradient  of  the  rate  of  buoy¬ 
ancy  generation  and  the  vertical  gradient  of  £  times  the 
tangential  component  of  frictional  stress.  It  follows  that 
surface  friction  can  induce  radial  motion  in  a  balanced 
formulation  of  the  boundary  layer. 

The  SE  equation  can  be  simplified  by  using  potential 
radius  coordinates  in  which  the  radius  r  is  replaced  by 
the  potential  radius,  R ,  defined  by 

\fR2  =  rv  +  1/r2, 

the  right-hand  side  being  the  absolute  angular  momen¬ 
tum  (Schubert  and  Hack,  1983).  In  this  case,  surfaces  of 
absolute  angular  momentum  are  vertical  and  the  assump¬ 
tion  that  these  surfaces  are  coincident  with  the  moist 
isentropes  provides  an  elegant  way  to  formulate  the  zero- 
order  effects  of  moist  convection  (Emanuel,  1986,  1989, 
1995,  1997).  It  is  for  this  reason,  perhaps,  that  bal¬ 
anced  models  remain  popular.  Nevertheless,  the  question 
remains  whether  the  boundary  layer  in  such  models  is 
sufficiently  accurate?  We  explore  this  question  below  in 
the  context  of  a  simple  slab  formulation  for  the  boundary 
layer  and  approximations  thereto. 


3.  A  slab  model  for  the  boundary  layer 

We  review  first  the  slab  boundary-layer  model  described 
by  Smith  and  Vogl  (2008;  henceforth  SV08),  which 
provides  a  suitable  framework  to  examine  the  accuracy  of 
various  approximations  in  boundary-layer  formulations. 
For  simplicity  we  assume  first  the  boundary  layer  to  have 
uniform  depth,  h ,  and  constant  density  as  in  Emanuel 
(1986).  In  our  cylindrical  coordinate  system,  the  vertically 
integrated  equations  for  the  radial  momentum,  azimuthal 
momentum  and  mass  continuity  can  be  written  in  the 
following  form: 
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where  Ub  and  Vb  are  the  vertically-averaged  radial  and 
azimuthal  components  of  wind  speed  in  the  boundary 
layer,  vg  (r)  and  Wh  are  the  tangential  wind  speed  and 
vertical  velocity  at  the  top  of  the  boundary  layer,  Cd  is 
the  surface  drag  coefficient  and  Wh-  =  ( Wh  —  \  wh\)/2. 
The  terms  involving  wsc  represent  turbulent  fluxes  at 
the  top  of  the  boundary  layer  (arising  from  rainbands, 
shallow  convection,  or  smaller-scale  turbulent  structures), 
but  for  simplicity  we  do  not  consider  them  here,  setting 
wsc  equal  to  zero.  Consistent  with  the  slab  formulation, 
the  quantities  Ub  and  Vb  are  assumed  to  be  independent 
of  depth.  Note  that  Wh-  is  non-zero  only  when  Wh  <  0, 
in  which  case  it  is  equal  to  Wh .  Thus  the  terms  involving 
Wh-  represent  the  transport  of  properties  from  above  the 
boundary  layer  that  may  be  different  from  those  inside 
the  boundary  layer.  As  in  SV08,  we  take  Cd  =  Cdo  + 
Cdi  |ub|,  where  Cdo  =  0.7  x  10-3  and  Cdi  =  6.5  x  10-5 
for  wind  speeds  less  than  20 ms-1  and  Cd  =  2.0  x  10-3, 
a  constant,  for  larger  wind  speeds.  These  values  are  based 
on  our  interpretation  of  Figure  5  from  Black  et  al.  (2007). 

Substitution  of  (11)  into  (9)  gives  an  expression  for  Wh'. 


Wh  = 


b+/(fg—  fb)  +  —  («b  +  vb)l^u^\ 


Ub 

r 


(12) 


where  a  is  zero  if  the  expression  in  square  brackets  is 
positive  and  unity  if  it  is  negative.  With  this  expression 
for  Wh ,  (9),  (10)  and  (12)  form  a  system  of  ordinary 
differential  equations  that  may  be  integrated  radially 
inwards  from  some  large  radius  R  to  determine  Ub  and 
Vb  as  functions  of  r,  given  values  of  these  quantities  at 
r  =  R  as  well  as  the  radial  profile  vg(r). 

In  the  following  subsection  we  examine  several 
approximations  to  the  foregoing  equations. 


3.1.  Scale  analysis 


At  this  stage  it  is  instructive  to  carry  out  a  scale 
analysis  of  (9)— (1 1).  (Whereas  some  of  the  criteria 
to  be  developed  are  implicit  in  the  modified  Oseen 
approximation  method  developed  by  Carrier  (1971)  for 
swirling  flow  boundary  layers  with  vertical  structure, 
the  scale  analysis  and  interpretations  herein  are  believed 
new.)  For  simplicity  we  consider  the  case  where  wsc  =  0. 
First  we  set  ub  =  vg  +  and  rewrite  the  two  momentum 
equations  in  the  form  given  in  Table  I.  There  we  have 
defined 


>ag 


dVa  Va 

=  -A  +  —  +  / 

dr  r 


to  be  the  absolute  vorticity  and  £g  =  2 vg/r  +  /  to  be 
twice  the  absolute  angular  velocity  of  the  gradient  wind 
profile  above  the  boundary  layer.  We  introduce  scales 
U,  V,  W  for  ( Ub ,  vg,  Wh~),  R,  H  for  the  radius  r  and 
boundary-layer  depth  h,  V '  for  v'h ,  A  for  £ag,  and 
S  for  £g.  The  two  terms  in  the  continuity  equation 
(11)  have  scales  U/R  and  W/H ,  respectively,  and 


Copyright  ©  2008  Royal  Meteorological  Society 


Q.  J.  R.  Meteorol.  Soc.  134:  1385-1395  (2008) 
DOI:  10.1002/qj 


1388 


R.  K.  SMITH  AND  M.  T.  MONTGOMERY 


Table  I.  Scaling  of  the  terms 

in  Equations  (9)  and  (10); 

see  text  for  details. 
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ub  —  = 
dr 

Wh- 
7  U  b 

h 

K2 
+ -!*- 

r 

+?gub 

Cd  ,  . 

—  Wb\ub 
h 

(9a) 

U2 

~R 

U2 

~R 

v'2 

~R 

ZV' 

> 

(9b) 

Su 

Su 

R°*v' 

1 

(9c) 

Terms  in  (10) 

dv' 

Wh-  , 

h Vb 

UbV b 

r 

£  ag^b 

-f/|Vb|(Vg  +  Vb) 

(10a) 

V' 

u— 

R 

V' 

w— 

H 

UV' 

R 

A  U 

—  V(V  +  V') 

H 

(10b) 

Rov > 

Rov' 

Rov' 

1 

V  V' 

Fvu  Fv~u 

(10c) 

Rov>  =  V'/(RA), 

Ro*v,  =  V'/(R 

3),  Su  = 

Ro*v,{U/V')2.. 

>  Fu  =  Cd 

V/(H  S)  and  Fv  =  CD 

V  /{HA). 

since  these  terms  have  equal  magnitudes,  it  follows 
that  W/H~U/R.  We  use  this  result  in  scaling  the 
two  terms  involving  Wh-  in  (9)  and  (10);  see  (9a)  and 
(10a)  in  the  table.  The  scales  of  each  term  in  these 
two  equations  are  given  in  the  first  row  below  each 
equation:  (9b)  and  (10b).  These  scales  are  made  non- 
dimensional  in  the  second  row  below  each  equation:  rows 
(9c)  and  (10c).  Five  non-dimensional  quantities  arise  in 
this  analysis: 

•  Rovr  =  V' /(R A),  a  local  Rossby  number  in  the 
tangential  momentum  equation  based  on  the  depar¬ 
ture  of  the  tangential  wind  from  the  gradient  wind 
(V7)  and  the  local  absolute  vertical  vorticity  of  the 
gradient  wind  above  the  boundary  layer  (A); 

•  Ro*,  =  V'/(RE),  a  local  Rossby  number  in  the 
radial  momentum  equation,  again  based  on  V',  but 
as  well  on  twice  the  absolute  rotation  rate  (based 
on  the  gradient  wind)  instead  of  A; 

•  Su  =  U2/(REV')  =  Ro*y,(U/V ')2,  a  local  Rossby 
number  multiplied  by  the  square  of  the  ratio  of 
agradient  winds,  ( U/V')2 ; 

.  Fu  =  CdV/(H 2)  and  Fv  =  CDV/(H A),  non- 
dimensional  forms  of  the  friction  terms.  These  are 
bulk-aerodynamic  Ekman  numbers  depending  on 
either  twice  the  absolute  rotation  rate  or  absolute 
vorticity,  respectively. 

Note  that  (U/V')2,  which  arises  in  the  expression  for 
S if ,  is  not  necessarily  small  and  might  be  larger  than 
unity. 

3.2.  Weak  friction  approximation 

Since  the  reduction  of  the  tangential  wind  component 
in  the  boundary  layer  increases  with  the  strength  of  the 
frictional  force,  it  is  of  interest  to  examine  first  the  case  of 
weak  friction,  for  which  \v^\  vg.  This  case  corresponds 
to  the  scaling  in  which  V' /V  1.  It  follows  immediately 
from  row  (10c)  in  Table  I  that 

U 

-  ~  Fv.  (13) 


Equation  (9a)  shows  that  even  in  the  case  of  weak  fric¬ 
tion,  the  net  radial  acceleration  —  §gu£  will  be  unop¬ 
posed  unless  it  is  balanced  by  the  radial  component  of 
the  frictional  force.  Such  a  balance  would  imply  that 
V'/V  ~  FU(U/V),  or 

V' 

y~FuFv.  (14) 

Clearly,  for  the  approximation  to  be  self-consistent,  the 
quantity  on  the  right-hand  side  of  (14)  must  be  much  less 
than  unity.  The  advective  part  of  the  radial  acceleration  in 
(9a)  will  be  negligible  locally  only  if  1,  requiring 
that 

U2  _  1  v'2 

V2  <<C  ~Rof, V2  (15) 

where  V'/V  is  obtained  from  (14).  Using  (13)  to  estimate 
U/V ,  it  follows  readily  that  (15)  is  equivalent  to  the 
condition  that  the  local  Rossby  number 

Roa  *«!.  (16) 

A  R 

As  shown  below,  this  condition  is  generally  not  satisfied 
in  the  high  wind  region  of  the  vortex  core  and  is  a  valid 
approximation  only  at  large  radii.  The  stringency  of  the 
condition  (15),  or  equivalently  (16),  will  become  apparent 
in  section  4.  The  second  term  on  the  right-hand  side  of 
(9b)  can  be  neglected  if  Ro*,  1,  i.e.  if  V'/V  Roe , 
which  we  define  as  V/(3^).  In  the  high  wind  speed 
region  of  the  vortex,  this  local  Rossby  number  is  about 
1/2.  Thus  in  this  region,  the  neglect  of  the  perturbation 
centripetal  acceleration  is  justified  if  V'/V  1/2. 

With  the  foregoing  scaling  assumptions,  (10a)  approx¬ 
imates  to  give 

ub  =  ~[ivg  ,  (17) 


where 
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Note  that  the  total  wind  in  the  boundary  layer,  |Vb|,  can 
be  approximated  as  vg  to  lowest  order  in  U/V  and  V'/V. 
At  this  level  of  approximation,  (9)  gives 

=  vb  -  vg  =  -v2vg ,  (19) 


where 


v  = 


Cdv  g 
hi 


(20) 


and  I2  =  £g£ag  is  the  inertial  stability  parameter.  Clearly 
the  integrity  of  the  approximations  leading  to  (17)  and 
(19)  depends  on  both  the  parameters  /x  and  v2  being 
small  compared  with  unity,  consistent  with  (15).  The  size 
of  these  parameters  may  be  readily  checked  for  a  given 
radial  profile  of  the  gradient  wind,  vg. 


3.3.  Semi-linear  approximation 

A  common  simplification  of  the  full  boundary-layer 
equations  is  to  neglect  the  terms  involving  the  vertical 
advection  through  the  top  of  the  layer  and  to  linearize 
the  acceleration  terms  (Eliassen  and  Lystad,  1977;  Kepert, 
2001).  If  we  make  such  approximations  in  (9)  and  (10), 
again  treating  the  magnitudes  of  u\>  and  as  small 
compared  with  vg9  but  not  linearizing  the  friction  term, 
the  equations  become 


Hvg-vb)  =  -^(u2  +  v2)V2ub,  (21) 

?  aMb  =  -  ^  (m^  +  vb  )  V2  ub  •  (22) 

n 

Again,  these  equations  may  be  solved  locally  for  u\>  and 
i>b  in  terms  of  vg,  £  and  f  a.  Dividing  (21)  by  (22)  gives 


(vg  ~  Vb)  2ub 

U  b  Vb 


(23) 


where  a2  =  f  ag/£g.  It  is  convenient  to  write  (wb,  Vb)  = 
vg (au,  v)  so  that  after  a  little  rearrangement,  (23)  gives 
the  equation  for  a  circle  in  the  (u,  v)  plane: 


u2  + 


1 

4' 


(24) 


o  O  o  o 

LD  o  LH  O  O 


U 


Figure  1.  The  circle  corresponding  with  Equation  (24)  and  the  curves 
given  by  Equation  (26)  for  typical  values  of  the  parameters  o  and 
v.  These  parameters  are  derived  from  the  radial  profile  of  vg  used  to 
construct  Figure  2.  The  straight  lines  are  those  given  by  Equation  (27) 
for  the  same  parameters.  The  numbers  marking  the  curves  are  the 
radii  (km)  at  which  the  parameters  a  and  v  are  calculated.  The 
darker  (blue)  curves  correspond  to  radii  <  250  km  and  the  lighter  (red) 
curves  for  radii  >  300  km.  This  figure  is  available  in  colour  online  at 
www .  interscience .  wiley .  com/j  oumal/ qj 


Families  of  curves  given  by  (24)  and  (26)  for  typical 
values  of  the  parameters  a  and  v  are  shown  in  Figure  1 
(see  section  4).  The  intersections  of  these  curves  in 
the  region  u  <  0,  v  >  0  (assuming  a  cyclonic  vortex  for 
vg)  give  the  required  solutions  for  ( u ,  v)  and  hence  for 
(u b,  fb).  These  solutions  may  be  readily  obtained  by  a 
simple  Newton-Rapheson  algorithm  as  explained  in  the 
appendix  to  SV08.  The  fact  that  all  solutions  for  (u,v) 
lie  on  a  circle  of  radius  one  half  implies  a  bound  on 
the  magnitude  of  the  ratios  cru  =  u\>/vg  and  v  =  Vb/vg, 
which  may  not  exceed  all  and  unity,  respectively. 
Therefore  the  occurrence  of  supergradient  winds  is  ruled 
out  in  these  linear  depth-averaged  formulations,  but  u\> 
could  exceed  vg  if  a  is  sufficiently  large. 

3.4.  Geostrophic  approximation 


Squaring  (22)  and  re-arranging  gives  a  second  equation 
relating  u  and  v: 

ti2v4  +  /jL2o2u2v2  —  o2u2  =  0  ,  (25) 

which  may  be  solved  v2  to  give 

i’=+a({i“4+^|  ~r2)  -  <26) 

noting  that  o2 n2  —  v2 . 


If  in  addition  to  the  semi-linear  approximation,  we  make 
the  small  Rossby  number  assumption,  £g  and  fag  may 
both  be  approximated  by  /  and  (21)  and  (22)  become 
the  slab  equivalent  of  the  Ekman  layer  equations  with 
a  quadratic  drag  law.  In  this  case  ar  =  1  and  (wb,  Vb)  = 
vg (u,  v).  Then  the  circles  given  by  (24)  remain  circles 
in  the  (u b,  Vb)  plane.  Moreover,  v  in  (26)  reduces  to 
vg  =  CDvg/(hf). 

3.5.  Linear  approximation 

Despite  the  ability  to  obtain  a  solution  to  the  approximate 
equations  with  the  quadratic  friction  terms  intact,  one 
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might  argue  that  a  more  consistent  approximation  of 
the  equations  would  include  also  a  linearization  of  the 
friction  terms,  replacing  ( u \  +  v by  vg,  but  retaining 
vb  on  the  right-hand  side  of  (22).  In  this  case,  (24)  does 
not  change,  but  (26)  is  replaced  simply  by 

v  =  --u.  (27) 


The  straight  lines  represented  by  (27)  are  shown  also 
in  Figure  1.  The  intersection  of  these  lines  with  the 
corresponding  circle  gives  the  solution  for  u  and  v  when 
the  quadratic  terms  in  (21)  and  (22)  are  linearized.  Note 
that  there  are  moderate  differences  between  the  solutions 
with  and  without  the  linearization  of  the  friction  terms 
for  most  values  of  the  parameters.  Now  (21)  and  (22) 
reduce  to 

Cd^ct 

$(vg  -  Vb)  = - b  ,  (28) 

h 

CDV* 

faMb  = - -*Wb.  (29) 

n 

These  equations  can  be  written  as 

vb  =  <72/JLU\)  ,  (30) 

and  ub  =  —fi(vg  +  vrh)  •  (31) 


They  have  the  solution 


and 

giving 


u  b  =  - 


1  +  v2’ 


vb 


2 

V  Va 
1  +  V2’ 


vb  = 


1  +  V2 


(32) 

(33) 

(34) 


According  to  the  scaling  analysis  in  subsection  3.2,  the 
approximation  is  self-consistent  only  if  /x2  v2 /Ro& 

1  /Ror.  Recall  that  Ro^  is  of  order  one-half  in  the 
vortex  core  region,  but  becomes  small  at  large  radius.  If 
/x  ^  v,  Su  ~  1  and  the  inertial  (advective)  terms  become 
important  in  the  radial  momentum  equation. 


3.7.  Emanuel’s  balanced  approximation 

Emanuel  (1986)  made  the  additional  assumption  that  the 
radial  wind  component  can  be  neglected  when  calculating 
the  total  wind  speed  in  the  friction  term.  In  this  case  the 
radial  wind  speed  is  given  by 


ub  =  — / uvg ,  (37) 

which  is  similar  to  the  dimensional  form  of  (27)  except 
that  vb  is  approximated  by  vg.  Moreover,  it  is  the 
same  as  (17),  which  is  a  feature  of  the  weak  friction 
approximation. 


4.  Some  solutions 

In  this  section  we  assess  the  accuracy  of  the  linear  and 
balanced  approximations  made  in  subsections  3.3,  3.4, 
3.5,  3.6  and  3.7  by  comparing  the  solutions  of  these 
equation  sets  with  the  corresponding  solution  of  the  full 
equations  (9),  (10)  and  (12),  which  we  refer  to  as  the 
benchmark  calculation  in  the  various  cases  discussed.  The 
profile  of  vg  is  the  same  as  that  used  in  SV08,  i.e. 

vg  =  v\s  exp  (— a\s)  +  V2 s  exp  (—012s), 

where  s  =r/rm,  and  rm  the  radius  at  which  the  tangential 
wind  speed  is  a  maximum  and  equal  to  vm.  In  the  calcula¬ 
tions  here,  vm  and  rm  are  taken  to  be  40  m  s-1  and  40  km, 
respectively,  and  the  other  parameters  are  chosen  to 
select  the  width  of  the  profile:  specifically  a\  =  1.4118, 
g?2  —  0.3,  v\  =  103.34ms-1  and  U2  =  20.0ms-1.  The 
Coriolis  parameter  is  taken  to  be  /  =  5  x  10-5  s-1.  The 
geostrophic  solution  is  used  to  initialize  the  benchmark 
calculation  at  a  radius  of  500  km,  where  the  local  Rossby 
number  vg{R)/{fR)  =  0.24  is  small  enough  for  this  solu¬ 
tion  to  be  reasonably  accurate  (here  vg(R)  =  5.88ms-1). 

4.1.  Constant-depth  boundary  layer 


3.6.  Balanced  approximation 

If  we  make  the  strict  gradient  wind  balance  approxima¬ 
tion  as  well  as  the  linear  approximation,  i.e.  if  we  assume 
that  the  tangential  wind  speed  in  the  boundary  layer  is 
equal  to  that  above  (ub  =  vg),  then  (21)  is  bypassed  and 
in  the  semi-linear  case  (22)  gives  directly  an  expression 
for  the  radial  wind  speed,  i.e. 

?a«b  =  --^(Mb  +  Ug)1/2ug'  (35) 

Squaring  this  equation  and  solving  for  ub  gives 


^b  —  1  • 

(1  -/X2)2 

The  balanced  approximation  cannot  be  rigorously  justi¬ 
fied  as  a  rational  approximation  using  the  foregoing  scale 
analysis  since  it  involves  setting  =  0  in  (9a).  Then 
the  value  of  ub  determined  by  (10a)  does  not  necessarily 
satisfy  the  remaining  terms  in  (9a). 


The  results  of  the  various  calculations  are  summarized 
in  Figure  2  for  a  boundary  layer  with  a  constant  depth 
of  800  m.  The  figure  shows  radial  profiles  of  (inward) 
radial  and  tangential  wind  components  in  the  boundary 
layer  and  the  vertical  velocity  component  at  the  top  of 
the  boundary  layer.  The  latter  is  calculated  analytically 
from  (12)  for  the  benchmark  calculation  and  from  the 
continuity  equation  (11)  for  the  approximations  to  it. 
Figure  2  shows  also  the  profile  of  tangential  wind 
assumed  at  the  top  of  the  boundary  layer  (vg)  and  the 
radial  variation  of  the  parameters:  \i  in  (17),  (18),  (32), 
(36)  and  (37);  v  in  (19)  and  (20);  o  in  (26);  and  v2  in 
(19),  (26)  and  (32)-(34). 

The  behaviour  of  all  solutions  shows  an  initial  increase 
in  r>b  and  \ub\  with  decreasing  radius  until  certain  radii, 
which  are  different  for  r»b  and  \ub\.  In  the  benchmark 
solution,  the  tangential  wind  speed  becomes  supergradient 
(i?b  >  Vg)  in  the  inner  core  region  (r  <  69  km)  and 
because  of  this,  the  radial  flow  suffers  a  rapid  deceleration 
(SV08)  and  declines  to  zero  at  a  radius  of  about  50  km.  At 
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r  (km)  r  (km) 


Figure  2.  Radial  profiles  of  (a)  tangential  and  (b)  radial  wind  components  in  the  slab  boundary-layer  calculation  described  in  section  3  and  in 
the  various  approximations  to  it.  The  solid  continuous  (black)  line  is  the  assumed  tangential  wind  profile  at  the  top  of  the  boundary  layer,  vg, 
and  those  that  terminate  near  a  radius  of  50  km  (red)  are  the  unapproximated  boundary-layer  solution.  Units  are  ms-1.  For  plotting  convenience 
the  sign  of  u b  has  been  reversed.  ‘SL’  represents  the  semi-linear  solution,  ‘L’  the  linear  solution,  ‘g’  the  geostrophic  solution,  ‘b’  the  balanced 
solution  and  ‘e’  Emanuel’s  balanced  solution,  (c)  The  corresponding  profiles  of  vertical  motion  at  the  top  of  the  boundary  layer,  (d)  The  radial 
variation  of  the  parameters  cr,  /jl ,  v,  v2  and  Roa  (labelled  ‘sigma’,  ‘mu’,  ‘nu’,  ‘nu2’,  and  ‘Ro’,  respectively)  for  the  assumed  tangential  wind 
profile.  This  figure  is  available  in  colour  online  at  www.interscience.wiley.com/joumal/qj 


this  radius  the  equations  become  singular  and  the  solution 
breaks  down.  Down  to  this  radius,  the  tangential  wind 
speed  is  still  increasing  inwards.  In  contrast,  in  all  the 
approximate  solutions,  the  tangential  wind  speed  remains 
subgradient,  attaining  a  maximum  and  then  decreasing. 
None  of  these  solutions  become  singular.  The  tangential 
wind  speed  in  the  semi-linear,  linear  and  geostrophic 
solutions  is  close  to  that  in  the  benchmark  calculation 
at  radii  larger  than  about  300  km,  but  the  geostrophic 
solution  shows  a  very  large  deviation  from  it  at  smaller 
radii.  In  contrast,  the  radial  wind  in  the  geostrophic 
calculation  follows  the  benchmark  rather  closely,  except 
at  radii  less  than  100  km,  whereas  that  in  the  semi-linear, 
linear  and  Emanuel  balanced  approximations  shows  a 
significant  deviation  from  the  benchmark  at  most  radii, 
being  a  considerable  overestimate  at  radii  between  about 
200  and  400  km  and  a  moderate  underestimate  at  radii 
around  100  km. 

The  balance  approximation  in  subsection  3.6  shows 
an  enormous  deviation  from  other  solutions,  a  feature 
that  is  attributed  to  the  fact  that,  at  least  for  the  profile 
chosen  for  vg,  the  parameter  /x  in  (36)  approaches 
unity  so  that  the  denominator  in  (36)  becomes  relatively 

Copyright  ©  2008  Royal  Meteorological  Society 


large  (Figure  2(d)).  Emanuel’s  balanced  approximation 
is  superior  in  this  respect  and  the  radial  wind  remains 
close  to  the  linear  solution  at  all  radii,  although  like  the 
latter,  it  shows  significant  departures  from  the  benchmark 
calculation.  There  is  very  little  difference  between  the 
radial  and  tangential  wind  components  in  the  semi-linear 
and  linear  calculations. 

There  are  considerable  differences  in  the  vertical  veloc¬ 
ity  between  the  various  approximations  and  the  bench¬ 
mark  calculation.  The  benchmark  calculation  breaks 
down  at  a  radius  of  about  50  km  where  the  radial  wind 
speed  tends  to  zero.  This  behaviour  is  not  replicated  by 
any  of  the  other  solutions.  At  slightly  larger  radii,  the 
vertical  velocity  in  the  benchmark  calculation  increases 
rapidly  with  decreasing  radius,  but  where  there  are  large 
radial  gradients,  the  approximations  on  which  boundary- 
layer  theory  is  based  will  no  longer  be  valid.  Again, 
the  balanced  solution  shows  a  large  deviation  compared 
with  all  other  solutions  between  about  200  and  400  km 
radius.  The  geostrophic  solution  predicts  large  vertical 
velocities  at  small  radii,  while  the  semi-linear,  linear  and 
Emanuel’s  balanced  solutions  show  a  maximum  vertical 
velocity  of  about  35cms-1  at  a  radius  of  about  60  km. 
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Figure  3.  Calculations  with  variable  boundary-layer  depth:  radial  profiles  of  (a)  assumed  boundary-layer  depth,  (b)  tangential,  and  (c)  radial  wind 
components  in  the  slab  boundary-layer  calculation  described  in  section  3  and  in  the  various  approximations  to  it.  (d)  shows  the  corresponding 
profiles  of  vertical  motion  at  the  top  of  the  boundary  layer.  Lines,  labels  and  units  are  as  Figure  2.  This  figure  is  available  in  colour  online  at 

www .  interscience,  wiley .  com/j  oumal/ qj 


There  are  significant  differences  also  between  the  radii 
at  which  the  vertical  velocity  changes  sign,  this  being 
about  170  km  in  the  benchmark  calculation,  235  km  in  the 
geostrophic  calculation,  300km  in  Emanuel’s  balanced 
calculation  and  310  km  in  the  balanced  calculation  of 
subsection  3.6.  The  scale  analysis  suggests  that  the  poor¬ 
ness  of  the  various  approximate  solutions  may  be  related 
to  the  assumption  that  Roa  1  and  /x2  v2 /Ro^ 
1/Roe  (section  3.5).  Figure  2(d)  shows  that,  at  least 
for  the  vg  profile  chosen,  these  conditions  are  far  from 
satisfied. 


5.  Solutions  for  radially  varying  layer  depths 

A  scale  analysis  of  the  full  boundary-layer  equations, 
without  vertical  averaging,  indicates  that  the  boundary- 
layer  depth  is  inversely  proportional  to  the  square  root 
of  the  inertial  stability,  /,  at  the  top  of  the  boundary 
layer.  This  scaling  would  imply  a  significant  reduction 
in  depth  between  the  starting  radius  and  the  core  region 
(e.g.  Figure  3(a)).  Such  a  variation  is  supported  by  the 
linear  solution  to  the  full  boundary-layer  equations  that 
accounts  for  vertical  structure  of  the  layer  (Eliassen  and 
Fystad,  1977;  Kepert,  2001)  as  well  as  full  numerical 
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solutions  (Kepert  and  Wang,  2001;  Montgomery  et  al ., 
2001).  While  it  is  not  possible  to  determine  the  radial 
variation  of  h  in  the  slab  model,  it  is  straightforward 
to  modify  (9)-(12)  to  allow  for  a  prescribed  variation 
h{r)  (see  SV08,  Appendix  A).  To  assess  the  effect  of 
a  decrease  in  h  with  declining  radius,  we  carried  out 
calculations  in  which  h{r)  =  h(R)^/(Ig/ /),  where  h{R) 
is  the  boundary  layer  depth  at  r  =  R  and  7g  is  the 
value  of  I  at  this  radius.  The  solutions  for  h{R)  =  800  m 
are  shown  in  Figure  3.  In  the  benchmark  calculation 
the  tangential  wind  speeds  in  the  boundary  layer  are 
decreased,  especially  inside  a  region  of  about  200  km  and 
the  peak  winds  are  significantly  lower  in  magnitude  than 
vm.  In  contrast,  the  peak  radial  winds  are  significantly 
larger  than  in  the  constant-depth  calculations  and  the 
maximum  occurs  at  markedly  smaller  radius.  When  the 
boundary-layer  depth  decreases  with  decreasing  radius, 
the  maximum  vertical  velocity  at  the  top  of  the  layer 
is  reduced  considerably  from  that  in  the  constant-depth 
calculations  and  is  more  in  line  with  that  in  previous 
calculations  (e.g.  Figure  3  of  Kepert  and  Wang,  2001). 
The  variable-depth  calculations  still  show  supergradient 
wind  speeds,  but  now  well  inside  rm  in  a  region  where 
radial  gradients  are  probably  steep  enough  to  strain  the 
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Figure  4.  As  Figure  3,  but  for  a  narrower  vortex  profile.  This  figure  is  available  in  colour  online  at  www.interscience.wiley.com/joumal/qj 


assumptions  of  boundary-layer  theory.  We  have  purposely 
omitted  to  show  the  balanced  solution  of  section  3.6  in 
Figure  3  as  the  parameter  /x2  in  (36)  has  values  larger 
than  unity  over  a  large  range  of  intermediate  radii.  At 
these  radii,  the  balanced  solution  does  not  exist  for  the 
particular  vg  profile  used. 

The  semi-linear  and  linear  approximations  predict  the 
tangential  wind  speed  rather  well,  almost  up  to  the  radius 
at  which  the  tangential  wind  component  becomes  super¬ 
gradient  in  the  benchmark  calculation.  The  geostrophic 
approximation  remains  accurate  inwards  to  a  radius  of 
about  300  km,  but  greatly  underestimates  the  tangential 
wind  component  inside  this  radius.  The  predictions  of  the 
radial  wind  component  by  the  approximate  theories  are  all 
poor  in  comparison  to  that  in  the  benchmark  calculation 
except,  surprisingly,  that  of  the  geostrophic  approxima¬ 
tion  at  radii  larger  than  about  200  km.  Inside  this  radius, 
however,  the  geostrophic  approximation  is  poor  also.  The 
Emanuel  balanced  approximation  and  hence  the  weak 
friction  approximation  give  by  far  the  worst  prediction 
of  the  radial  component  with  a  maximum  exceeding  that 
of  the  benchmark  calculation  by  10-20 ms-1  at  radii 
between  about  70  and  300  km.  Note  that  in  these  approx¬ 
imations,  the  radial  wind  is  just  minus  the  ratio  of  the 
tangential  wind  above  the  boundary  layer  multiplied  by 
the  profile-dependent  parameter  \i.  All  the  approximate 
theories  overestimate  the  subsidence  into  the  boundary 
layer  at  outer  radii  (>  320  km)  and  show  ascent  occurring 
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inside  this  radius  compared  with  a  radius  of  160  km  in 
the  benchmark.  Moreover  the  ascent  is  much  larger  than 
in  the  benchmark  except  within  a  radius  of  about  30  km. 

Figure  4  shows  plots  analogous  to  those  in  Figure  3,  but 
for  a  narrower  vortex  profile  that  has  the  same  values  of 
rm  and  vm.  (The  functional  form  is  the  same  as  before,  but 
v\  =  90.06ms-1,  V2  =  36ms-1,  a\  =  1.601,  oii  =  0.5.) 
The  radius  of  gale-force  winds  (17ms-1)  is  317km 
compared  with  381  km  in  Figure  3.  For  this  profile,  the 
inertial  stability  parameter,  /,  has  a  local  minimum  at  a 
radius  of  about  240  km,  which  implies  a  local  maximum 
in  the  boundary-layer  depth  at  this  radius  (Figure  4(a)). 
While  this  maximum  is  a  natural  consequence  of  the 
scaling,  it  is  not  known  how  realistic  it  is  because  of  the 
lack  of  observational  data  on  the  depth  of  the  inflow  layer 
at  these  radii  in  a  tropical  cyclone.  Because  the  radial 
gradients  are  larger,  one  would  expect  a  narrow  vortex 
profile  to  strain  the  balanced  approximations  more  than 
a  broader  one  and  this  expectation  is  confirmed  by  the 
solutions.  Now  the  semi-linear  balanced  approximation 
is  the  worst  of  all  of  the  approximations  in  capturing 
the  tangential  wind  component  in  the  boundary  layer  at 
radii  between  60  and  300  km.  In  contrast,  the  semi-linear 
and  linear  calculations  are  reasonably  accurate  until  near 
the  radius  at  which  the  tangential  wind  component  in  the 
benchmark  calculation  becomes  supergradient.  However 
the  radial  wind  component  is  greatly  overestimated  by  all 
the  approximations  at  a  radii  less  than  about  320  km,  with 
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Figure  5.  Contributions  to  the  radial  momentum  balance  for  radii  less  than  200  km  in  the  variable-depth  calculations  in  (a)  Figure  3,  and  (b)  the 
narrower  vortex  profile  in  Figure  4.  The  labelling  on  the  curves  refers  to  the  net  inward  force  (nif),  the  frictional  force  (fri),  and  the  vertical 
advection  of  momentum  into  the  boundary  layer  (wu).  The  residual  force  is  the  solid  (blue)  line.  This  figure  is  available  in  colour  online  at 

www .  interscience,  wiley .  com/j  oumal/ qj 


the  Emanuel  balanced  calculation  being  the  least  accurate. 
The  vertical  velocity  at  the  top  of  the  boundary  layer  is 
likewise  very  poor  compared  with  that  in  the  benchmark 
calculation  and  worse  than  in  the  broader  tangential  wind 
profile  in  Figures  2  and  3. 

The  reasons  for  the  large  deviations  of  the  radial  and 
vertical  winds  in  the  balance  solutions  from  those  in  the 
benchmark  calculation  may  be  inferred  from  Figure  5, 
which  shows  the  radial  force  fields  for  the  two  vortex 
profiles  in  the  calculations  with  variable  boundary-layer 
depth.  The  net  inward  driving  force,  the  difference 
between  the  inward  pressure  gradient  force  and  the 
outward-directed  centrifugal  and  Coriolis  forces,  exceeds 
the  frictional  force  by  a  significant  fraction  leaving  a  net 
inward  force  which  is  just  the  radial  acceleration.  This  net 
inward  force  is  appreciable  down  to  the  radius  just  inside 
rm,  when  it  begins  to  decrease  rapidly,  reversing  sign 
about  10  km  outside  the  radius  at  which  the  net  driving 
force  changes  sign  and  the  tangential  wind  component 
becomes  supergradient.  The  foregoing  remarks  apply  to 
both  panels  in  Figure  5  and,  as  expected,  the  radial 
acceleration  is  a  little  larger  for  the  narrower  vortex 
profile. 

It  is  worth  pointing  out  that  the  benchmark  calculation 
itself  is  not  entirely  realistic  in  these  varying-depth  calcu¬ 
lations.  For  one  thing,  the  150  m  depth  near  the  radius  of 
maximum  gradient  wind  speed  is  a  factor  of  two  to  three 
too  shallow  in  the  inner  core  region  when  compared  to 
observations  (e.g.  Bell  and  Montgomery,  2008;  Kepert, 
2006a, b)  and  for  another,  the  predicted  inflow  angles, 
tan-1(|w|/i;),  are  too  large,  on  the  order  of  45°,  com¬ 
pared  with  observed  maxima  of  typically  20-30°  (e.g. 
Kepert,  2006a, b).  These  two  shortcomings  are  presum¬ 
ably  related.  The  shallow  depth  is  a  direct  consequence 
of  the  simple  Elias  sen  and  Fystad  scaling  for  a  bound¬ 
ary  layer  with  continuous  vertical  variation  assuming  a 
constant  eddy  diffusivity  K  and  could  be  addressed  by 
allowing  K  to  increase  inwards.  Despite  these  shortcom¬ 
ings,  we  believe  that  the  results  are  interesting  and  worth 
reporting. 
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6.  Conclusions 

Our  calculations  for  the  chosen  tangential  wind  profiles 
vg  show  that  the  semi-linear  and  linear  solutions  as  well 
as  Emanuel’s  balanced  solution  (i.e.  the  weak  friction 
solution)  tend  to  overestimate  the  inflow  in  the  bound¬ 
ary  layer  at  large  radii  (>  160  km)  and  underestimate  it 
at  inner  radii.  In  some  cases  they  predict  the  maximum 
vertical  motion  to  occur  at  an  unrealistically  large  radius 
(>250  km)  and  they  are  unable  to  capture  the  formation  of 
supergradient  winds.  The  balanced  formulation  without  a 
linearization  of  the  friction  term  is  especially  unrealistic. 
The  inaccuracy  of  these  approximations  is  suggested  by 
a  rational  scale  analysis  of  the  slab  boundary-layer  equa¬ 
tions  since  the  conditions  under  which  the  approximations 
are  derived  are  not  satisfied,  at  least  for  the  parameters 
and  gradient  wind  profiles  studied  here.  Therefore,  we 
consider  such  formulations  to  be  inappropriate  for  repre¬ 
senting  the  boundary  layer  in  the  inner  core  region  of  a 
hurricane. 
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